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ABSTRACT 

Low- ionization transitions such as the Mgn A2796/2803 doublet trace cold gas in the vicinity of 
galaxies. It is not clear whether this gas is part of the interstellar medium of large proto-disks, part of 
dwarfs, or part of entrained material in supernova-driven outflows. Studies based on Mgn statistics, e.g. 
stacked images and clustering analysis, have invoked starburst-driven outflows where Mgn absorbers 
are tracing the denser and colder gas of the outflow. A consequence of the outflow scenario is that the 
strongest absorbers ought to be associated with starbursts. We use the near-IR integral field spectrograph 
SINFONI to test whether starbursts are found around z ~ 1 Mgn absorbers. For 67% (14 out of 21) 
of the absorbers with rest-frame equivalent width W r A2796 > 2 A, we do detect Ha in emission within 
±200 km s _1 of the predicted wavelength based on the Mgn redshift, and with impact parameter ranging 
from 0.2" to 6.7" from the QSO. The star formation rate (SFR) inferred from Ha is 1-20 M yr^ 1 , 
i.e. showing a level of star formation larger than in M82 by a factor of > 4 on average. Our flux limit 
(3-cr) is /jj a < 1.2 x 10~ 17 erg s _1 cm~ 2 , corresponding to a SFR of ~0.5 M yr _1 , much below past 
ground-based Ha surveys of absorption-selected galaxies. We find evidence (at > 95%) for a correlation 
between SFR and equivalent width, indicating a physical connection between starburst phenomena and 
gas seen in absorption. In the cases where we can extract the velocity field, the host galaxies of Mgn 
absorbers with W^ 2796 > 2 A reside in halos with mean mass < log M^(M Q ) >~ 11.2 in good agreement 
with clustering measurements. 

Subject headings: cosmology: observations — galaxies: high-redshift — galaxies: evolution 



1. INTRODUCTION 

In the era of large spectroscopic surveys such as 
SDSS (York et al. 2000), the Mgn A2796/2803 doublet 
seen in background QSO spectra is particularly useful for 
selecting thousands of absorbing galaxies (e.g. York et al. 
2006). Galaxies selected via their Mgn A2796/2803 dou- 
blet have been used as proxies for studying Hi selected 
galaxies at z < 1.5 (e.g. Rao et al. 2006). 

It is not clear whether this cold gas is part of the inter- 
stellar medium of large spirals (Wolfe et al. 1986), part of 
the halos of galaxies (Bahcall & Spitzer 1969) or part of 
dwarf galaxies (York et al. 1986) akin to the Magellanic 
Clouds. An alternative scenario for absorption-selected 
galaxies is that the gas seen in absorption is part of cold 
gas clumps in the host galaxy halos entrained in outflows 
produced by supernovae (SNe) (Nulsen et al. 1998; Schaye 
2001). Most likely all three scenarios play a role, but in 
what proportion for a given equivalent width or Hi column 
density? 

Using cross-correlation techniques, Bouche et al. (2006) 
statistically constrained the halo mass of ~2000 Mgn host 
galaxies, finding that the host halo mass (M^) decreases 
with increasing Mgn rest-equivalent width (IF/ 12796 ). 
Given that the equivalent width must be correlated with 



the line-of-sight velocity width, Av, as also verified obser- 
vationally by Ellison (2006), these results imply that Mh 
and Av are anti-correlated. This seems to show that the 
clouds responsible for the absorption are not virialized, 
otherwise, a M^—Av correlation would have been mea- 
sured. A natural explanation for this is that the cold gas 
originates from SN-driven outflows, in which case the ve- 
locity width Av may be related to the outflow kinematics 
or may be a measure of the mass outflow rate. In a few 
cases, evidence for super- winds is seen from the absorption 
profile (Bond et al. 2001; Ellison et al. 2003). Others, e.g. 
Prochter et al. (2006) and Martin (2006) arrived at the 
same conclusion that super-winds play a significant role in 
Mgn absorbers from very different perspectives. 

The most important implication from the starburst sce- 
nario is that the selection of galaxies via the Mgll absorp- 
tion signature may be equivalent to selecting starburst- 
producing super- winds. Equivalently, non-detections of 
the Ha signature of a starburst would enable one to rule 
out the starburst scenario, while the presence of Ha emis- 
sion is no un-ambiguous proof of starburst-driven winds. 

In this letter, we report the first results of our "SIN- 
FONI Mgn Program for Line Emitters" (SIMPLE), which 
is aimed at detecting Ha from the starburst galaxy using 



1 Based on observations made at the ESO telescopes under program ID 077.A-0576, 078.A-0600, 078.A-0718 and 079.A-0600. 

2 Max Planck Institut fur extraterrestrische Physik, Giessenbachstrasse, D-85748 Garching, Germany; NB: nbouche@mpe.mpg.de 

3 Centre for Astrophysics & Supercomputing, Swinburne University of Technology, Hawthorn, Victoria 3122, Australia 

4 European Southern Observatory, Karl-Schwarzschild-str 2, D-85748 Garching, Germany 



1 



2 



the integral field unit (IFU) SINFONI (Eiscnhauer et al. 
2003) available at the Very Large Telescope (VLT). The 
advantages of IFUs include the possibility of detecting the 
Mgll host at impact parameters smaller than the seeing 
disk, and of measuring the two-dimensional kinematics of 
galaxies. 

We present our sample selection in section 2, and results 
in section 3. Throughout, we use the h = 0.7, Dm = 0.3, 
Oa = 0.7 cosmology. 

2. SAMPLE SELECTION 

Our sample of Mgn-selected galaxies is selected from 
the SDSS/2QZ databases and from the compilation of 
Ryabinkov et al. (2003) with the only physical criterion 
jyA2796 > 2 A. We then ensured that the corresponding 
Ha emission line would fall inside the SINFONI wave- 
length range and away from sky OH emission lines. The 
W^ 2796 > 2 A criterion was used for the following two rea- 
sons. First, in the super-wind scenario, the strongest ab- 
sorbers (as measured by W r A2796 ) ought to have the largest 
star-formation rates, hence the largest Ha fluxes. Second, 
the W^ 2796 > 2 A criterion is known to select the hosts 
with the smallest impact parameters, p < 30 h~ x kpc (e.g. 
Steidel 1995; Bouche et al. 2006). At z ~ 1, this corre- 
sponds to ~ 4", which means that the host galaxy will fall 
within the field of view of SINFONI (8"). 

We have obtained SINFONI observations towards 21 
Mgll absorbing galaxies all with absorption redshift z a b s — 
1, i.e. corresponding the J-band filter. We focused on 
z <~ 1 systems because (1) this redshift range is close to 
the range used in Bouche et al. (2006), and (2) our flux 
sensitivity is enhanced in the J-band compared to that in 
the K-band. The latter point is due to the gain in lu- 
minosity distance compared to the loss in raw instrument 
throughput. For each QSO field we took 4x600s expo- 
sures, placing the QSO in the four quadrants of the field 
of view. The total exposure time for the host galaxy thus 
varies between 10 and 40 minutes, depending on its loca- 
tion with respect to the QSO. The observations were all 
taken in secing-limitcd mode with 0.125" pixels, in moder- 
ate (optical) seeing conditions (1"— 1.5"), yielding a point 
spread function (PSF) of 0.8"-1.2" in the' J-band. 

The data reduction is based on pairwise subtraction of 
the frames as in Forster Schreiber et al. (2006) using the 
MPE SINFONI pipeline (SPRED, Abuter et al. 2006, and 
references therein) and the OH sky line removal scheme 
of Davics (2007). The wavelength solution is obtained us- 
ing the Ar lamp frames, and the absolute wavelength is 
obtained directly from the OH sky emission lines. Flux 
calibration is performed using telluric standard stars. 

3. RESULTS 

From the SINFONI data-cubes, we detect an Ha emit- 
ter within ±200 km s _1 of the predicted location from the 
Mgll redshift in 14 of the 21 fields down to a flux limit of 
4xl0~ 18 erg s _1 cm~ 2 (1-cr) (see Table 1). In other words, 
we have a 67% detection rate. The undetected hosts are ei- 
ther fainter than our flux limit or outside the field of view. 
Fig. 1 shows the Ha flux maps for all of the detected host 
galaxies. In all frames, the QSO PSF (typically 0.8-1.2" 
FWHM) was subtracted from the Ha flux maps and is 
shown as contours. 



As seen from Fig. 1, the impact parameters range from 
0.2" to 4" (1 to 25 kpc) with one at ~ 6.7" (J0841+2339). 
The very strong absorber (W / r A2796 =4.5A) towards the 
QSO 2QZJ0226-28 is detected with an impact parameter 
of 0.2", well within the seeing disk. We can rule out the 
possibility of this emission line being from the QSO. In- 
deed, the emission line at 1.327^im corresponds to 4185 A 
in the rest-frame of the QSO and therefore does not cor- 
respond to any known atomic transition. Furthermore, 
the emission line is clearly spatially offset from the QSO 
centroid. This shows that it is possible to decouple the 
intervening galaxy and the QSO owing to the information 
contained in the IFU data cube. 

It is extremely unlikely that all the 14 emission lines 
shown in Fig. 1 be something else in relation to other in- 
tervening systems along the line of sight. For the strongest 
emitters (J0448,J0302,J0822,J0943), we do also detect 
[Nil] A6583 and [SII] AA6717,6731, therefore lifting any 
ambiquity. In the cases where only one line is detected, 
we can reject the other possibilities (On, Om, H/3) for the 
following reasons. [Om] A5007 is a doublet easily identi- 
fiable at the SINFONI resolution R ~ 2200. [On] A3727 
would not be resolved, but that would put the galaxy at 
a redshift higher than the QSO emission redshift. In most 
cases, the same argument applies for H/3. When not, we 
have checked the QSO spectrum for possible correspond- 
ing absorbers (Mgll, Fell) and found nothing in the QSO 
spectrum at the corresponding location. 

After more than two decades, only < 100 absorption- 
selected galaxies have been identified in emission, and 
mostly at z ~ 1 (e.g. Kulkarni et al. 2006, and refer- 
ences therein) using broad band imaging (e.g. Le Brun 
et al. 1997; Chen & Lanzetta 2003). Low-rcdshift (z < 1) 
Mgll absorbers have been identified more easily (Bergeron 
& Boisse 1991; Steidel et al. 1994; Le Brun et al. 1997) 
than low-z and high-z damped Lya absorbers (DLAs), 
which are all also Mgll absorbers (Rao et al. 2006). 
Our 67% detection rate of the host Ha emission is thus 
in contrast to most previous Ha studies of absorption- 
selected galaxies since our sample of strong Mgll ab- 
sorbers (W^ 2796 > 2 A) is dominated by DLA 'candi- 
dates' meeting the criterion of Rao et al. (2006) (see Ta- 
ble 1) and is due to several factors. If our careful avoid- 
ance of OH sky emission lines plays a significant role, it 
is the SINFONI throughput that enabled us to reach flux 
limits unavailable previously from the ground, reaching 
1.2 x 10~ 17 erg s" 1 cur 2 or a SFR of -0.5 M Q yr" x at 
z = 1 in less than 1 hr. For comparison, previous Ha sur- 
veys for z > 2 DLAs include the ground-based surveys 
of (e.g.) Mannucci et al. (1998) unveiling several can- 
didates down to 1.4 x 10~ 16 erg s _1 cm~ 2 (3cr), van der 
Werf et al. (2000) unveiling several candidates down to 
1.7 x 10~ 16 erg s" 1 cur 2 (3er), Bunker et al. (1999) show- 
ing no detection down to 7-16 x 10 _17 erg s _1 cm~ 2 (3cr), 
the fiST/NICMOS survey of the z = 0.656 DLA 
(W r A2796 =1.5 A) by Bouche et al. (2001) showing no detec- 
tions down to 3.7 x 10 _17 erg s _1 cm~ 2 (3er), and the NIC- 
MOS survey of a z = 1.89 DLA by Kulkarni et al. (2001) 
showing no detection down to 1.4x 10~ 17 crg s _1 cm~ 2 (3cr). 

The large number of Ha detections in this survey shows 
the relatively high star formation rates of the host galaxy 
given our shallow exposures. The averaged integrated Ha 
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luminosity for our sample is ijj a — 6.4 x 10 erg s 1 
(median 4.2 x 10 41 erg s _1 ), uncorrected for extinction. 
This is about > 4 times that of the starburst M82 whose 
raw L Ha is 1.5 x 10 41 erg s _1 (Lehnert et al. 1999). If 
we adopt a statistical correction of Ay = 0.8mag as de- 
rived for local star-forming and starburst galaxies of sim- 
ilar Ha luminosities (Kennicutt 1998; Buat et al. 2002), 
we infer an average (median) SFR of <~12 (7) M Q yr -1 
using the Kennicutt (1998) flux conversion for a Salpeter 
IMF, and -7 (4) M© yr" 1 using the Chabrier (2003) IMF 
(both from 0.1 M Q to 100 M ). Table 1 lists the derived 
properties of our sample. 

Fig. 2 shows the SFR as a function of W,. 796 . It appears 
that there is a correlation between the SFR and W^ 2796 , 
which is significant at more than 2-cr: the Pearson correla- 
tion coefficient is p — 0.59 with a P-valuc= 0.024 < 0.05, 
and a Spearman rank test gives p = 0.53 and P = 0.048. 
The corrrelation is very consistent with the analysis of Zi- 
betti et al. (2007) who found that the hosts of the strongest 
Mgn absorbers are bluer. If stronger starbursts produce 
larger outflows (whose signature is more clouds over a 
larger Av), such a correlation would naturally arise in the 
outflow scenario. We will attempt to model this in more 
detail in Bouche et al. (in preparation). 

The inferred SFRs are moderate, but whether or not the 
Ha emitters seen in Fig. 1 are producing winds depends on 
the SFR per unit area, E, not on the integrated SFR. For 
instance, the starbursts M82 (SFR— 2 M Q yr -1 , Lehn- 
ert et al. 1999) and NGC 1569 (SFR- 1 M Q yr" 1 , Mar- 
tin, Kobulnicky & Heckman 2002) have E greater than 
the 0.1 M Q yr -1 kpc~ 2 threshold for harboring SN driven 
outflows (Heckman 2003), and do show strong outflows. 
Already, the majority of our sample meet this criteria. 
Higher spatial resolution data and better size constraints 
will enable us to test whether strong Mgn selected galaxies 



always meet the outflow threshold and, in the cases where 
it does, to probe the outflow physical properties. 

As stated in section 1, Ha ought to be present and 
detected in the starburst-outflow scenario. Had we not 
detected Ha, this would have ruled out the outflow sce- 
nario. The high detection fraction thus prevents us from 
ruling out the outflow scenario. The relatively high av- 
erage SFR (compared to M82) is in fact very consistent 
with the prediction of the starburst-outflow scenario, and 
we view this result as supporting the outflow scenario, but 
other interpretations are possible. Note that our conclu- 
sion only applies for galaxies selected with the strongest 
Mgn absorption equivalent widths (W r A2796 > 2 A). Mgn 
absorption systems with W^ 2796 < 2 A are much more 
numerous and may select a more diverse class of galaxies 
(e.g. Bergeron & Boisse 1991; Steidel et al. 1994). 

Fig. 3 shows an example of the two dimensional veloc- 
ity field for one of the Mgn host galaxies. For the re- 
mainder of the sample, when we can extract the veloc- 
ity field, we find that the host galaxy of Mgn absorbers 
with VF r A2796 > 2 A reside in halos with masses rang- 
ing from log M/j(M ) — 9.9-11.9 with a logarithmic mean 
< log M/j(M ) >— 11.2, once corrected for inclination. 
This is 10 times smaller than the Milky Way, and in good 
agreement with the halo masses predicted from the cluster- 
ing analysis of Bouche Murphy Peroux (2004) and Bouche 
et al. (2006). Here we are using the maximum velocity 
Knax as a proxy for the circular velocity. In N. Bouche et 
al. (in preparation), we will present the kinematic results 
in detail. 
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SDSSJ014717. 76+125808. 8 EW=4.0 2QZJ022620. 4-285751 EW=4.5 2QZJ024824.4-31 0944 EW=2.5 2QZJ030249. 6-321 600 EW=2.2 J0427-1302 EW=2.0 




Fig. 1. — We show the 14 QSO fields (sorted by RA.) where the Mgn host galaxy was detected via its Ha emission using the IFU SINFONI. 
The contours show the QSO continuum (showing the field seeing). The arrow represents 1" (~ 8 kpc at z ~ 1) and the impact parameters 
range from 0.2 to 4". One galaxy is seen 0.2" away from the QSO (2QZJ0226-28), well inside the seeing disk (0.8") shown by the QSO 
continuum contours, showing the power of untangling spatially overlapping light when using the full 2-dimcnsional spectral information from 
IFU data. North is up, East is left, color-coding reflects the relative Ha fluxes of the sources. 
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Fig. 2. — The star-formation rate (SFR) from the Ha emission of Mgn hosts as a function of rest-frame equivalent width W™ 796 . SFR 
and W,- A2796 may be correlated, indicating a physical connection between the host starburst and the gas seen in absorption. The cut at 
w \2796 = 2 A is due to our sample selection. Our SFR flux limit (SNR=3) is shown as the dotted line. 



Fig. 3. — Comparison between the Ha emission map (left) and the two dimensional velocity field (right) for one of our Mgn host galaxies. 
The QSO continuum is shown as grey contours. The kinematics reveal that the kinematic major axis is perpendicular to the QSO-galaxy 
vector, i.e. the QSO is located along the minor axis, which is further evidence that the gas traced by Mgn is part of the halo gas and not 
the disk. The maximum rotation velocity is V ma x ~ 80 km s — 1 for this galaxy and is typical for our sample. Using V max as a proxy for the 
circular velocity, the inferred mean halo mass is < logMh(AfQ) >~ 11.2 for the sample. 
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1.68601.02094 3.608 


1 


1.1 


NA 


NA 


< 1.2E-17 


< 0.6 


2QZJ232330.4-292123 


1.547 0.811 2.0430 




1.0 


NA 


NA 


< 1.2E-17 


< 0.6 



Table 1 
Sample Properties. 



a a '1' indicates that the system meets the Rao et al. (2006) criteria for being a DLA. 

6 the 3-er upper limits for the non-detections are computed for an unresolved source spread over ~ 200 pixels, spread over 32 spatial pixels and 
spectral FWHM=6 pix=9A (i.e. R 1/2 = 0A"~ 3 kpc, and FWHM=250 km s" 1 ). The same dust-correction was applied to the SFR limits. 



